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ABSTRACT 

We use deep multi-band optical and near-infrared data for four general fields, GOODS-South, HDF 
North/South, and IRAC UDF in GOODS-North to investigate the evolution of the observed rest- 
frame U — V color of field galaxies as a function of the stellar mass evaluated by fitting the galaxy 
spectral models to the observed broad-band SEDs. In these four fields, the U — V color distributions 
of the galaxies at each mass and redshift interval are very similar. At 0.3 < z < 2.7, we found that 
more massive galaxies always tend to have a redder U — V color. High- and low-mass galaxies exhibit 
quite different color evolutions. As seen in our previous study in HDF-N, the color distribution of 
low-mass (M ste uar S 3 x 10 9 Mq) galaxies becomes significantly bluer with an increase in the redshift. 
This evolution of the average color can be explained by a constant star formation rate model with 
Zform ~ 4. On the other hand, the average color of high-mass galaxies (M ste u a r ^ 3x 10 10 M ) evolves 
more strongly at a high redshift. Such mass-dependent color distribution and its evolution indicate 
that galaxies with a larger stellar mass appear to have shorter star-formation timescales, and on an 
average, they form the larger fraction of their stars in the earlier epoch. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

Stellar mass is one of the most important quan- 
tities of galaxies that facilitate the understanding of 
their formation and evolution. The evolution of the 
number density of galaxies as a function of the stel- 
lar mass conveys the overall history of star forma- 
tion and the mass assembly process in the universe. 
Many authors have investigated such evolution of the 
stella r mass distribution function back to z > 1 
e.g.. iDickinson. Papovich. Ferguson. &: Budavari 



Bell et al.l 120031 iFontana et all 12004 IDrorv etal 
Drorv et al.ll2005fl 



2003, 



2004, 



On the other hand, the study of properties such as 
the spectral energy distribution (SED), colors, and line 
fluxes of galaxies as a function of the stellar mass at var- 
ious epochs helps us to understand star formation his- 
tories in additional detailed, such as, in what manner 
the star formation occurs in low- and high-mass galaxies 
and how it contributes to the growth of their stellar mass. 
In the case of l ocal galaxies d e tected in t he Sloan Digi- 
tal Sk y Survey. Illeavens et alJ l)2004|) and lJimenez et all 
( 2005) argued that most of the stars in the present-day 
massive galaxies have been formed at a high redshift, 
while star formation histories of low-mass galaxies have 
relatively longer timescal es. In more d i rect o bservations 
of high-redshift galaxies, ICowie et alJ ljl996fl originally 
found that the mass of the galaxies undergoing rapid 
star formation has been declin ing smoothly with a de- 
crease in the redshift. Recently. I Juneau etal] l|200l also 
investigated the stellar mass dependence of star forma- 
tion rates of galaxies at 0.8 < z < 2.0 using the Gemini 
Deep Deep Survey data. They found the trend to be 
similar to that predi cted in the studies of local galaxies 
ijHeavens et alJ2004j) . Similar results are obtained in the 

Electronic address: kajisawa@optik.mtk.nao.ac.jp 



studies on the ev olution of specific star formation rates 
back to z ~1-1.5 tBrinc hmann fc Ellisll200fl iBauer et alJ 
120051 IFeulner et aUl2005aj) . Very recently. iFeulner et alJ 
( 2005bJ) extended the analysis to z ~ 5, and found that 
a similar trend continues even at a higher redshift. 

The samples in these previous studies, however, are 
too shallow to reach a sufficiently low mass range be- 
low ~ 5 x 10 9 Mq, a characteristic "transition" mass 
(Kajisawa and Yamada 2005, and Section 4 below), or 
are biased toward objects with typically low M/L ratios, 
such as more actively star-forming bluer populations due 
to the lack of very deep near-infrared (NIR) data. The 
evolution of low-mass galaxies is, however, important to 
understand the overall picture of galaxy evolution not 
only because they may be the building blocks of more 
massive galaxies but also because they are more sensitive 
probes to understand the thermal history of the universe. 

In Kajisawa and Yamada (2005, hereafter referred to 
as KY05), we investigated the evolution of the rest- 
frame color and morphology distribution of the galaxies 
in HDF-N as a function of the stellar mass back to z ~ 3 
over a wide range of stellar masses down to 10 9 M Q . We 
found that the behavior of the rest-frame U — V color dis- 
tribution and its dependence on the stellar mass change 
at around ~5x 10 9 Mq. In addition, we discovered that 
the sequence of the color distribution of low-mass galax- 
ies systematically becomes bluer with the redshift, while 
the colors of high-mass galaxies do not show a signifi- 
cant redshift evolution. However, the area of the HDF-N 
data used in our analysis was only ~ 3.5 arcmin 2 , and it 
was inferred that the effect of field-to-field variation was 
large. 

Therefore, in this paper, we have extended the study 
area of KY05 to inclu de some other fields, GOODS-S 
ijGiavalisco et al.ll2004fl. a part of the IRAC Ultra Deep 
Field in GOODS-N HfDickinson et"ai]l2003bj) . and FIRES 
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fHDF-South. lLabbe et afll2003|) . Using the deep optical 
and NIR data for these fields, we investigate the rest- 
frame U — V color evolution of galaxies at 0.3 < z < 
2.7 as a function of the stellar mass. These data allow 
us to extend our previous analysis to a larger sample 
in order to study the evolution of galaxies with a high 
statistical accuracy. We used the Vega magnitude system 
throughout the paper. 

2. DATA & SOURCE DETECTION 

We used multi-band optical and NIR data in the four 
fields, GOODS-S, FIRES, HDF-N, and a part of IRAC 
UDF in GOODS-N. For HDF-N, we used the same data 
as those used in KY05 (HST/WFPC2, HST/NICMOS, 
Subaru/CISCO, UBVIJHK set). For GOODS-S, we used 
the publicly available VLT/ISAAC verl.5 data (JHK- 
bands) and HST/ACS verl.O data (bviz, Giavalisco et 
al. 2004). For FIRES, we used the publicly available 
HST/WFPC2 + VLT/ISAAC data (UBVIJHK, Labbe 
et al. 2003). 

For IRAC UDF in GOODS-N, we observed a field in 
the J and K bands using the Subaru Telescope equipped 
with CISCO (Motohara et al. 2002) on March 10, 2003. 
The field of view (FOV) used in the analysis was 1.74' 
x 1.74', and the total exposure time in the K and J 
bands was 14780 and 3880 s, respectively. The data were 
reduced using the IRAF software package, as described in 
KY05. The publicly available GOODS HST/ACS verl.O 
data (bviz) were also used in the field. 

For all the fields, we independently performed source 
detection in the iVban d images using the SExt ractor 
image analysis package l)Bertin fc Arnoutsl IT996). We 
adopted MAG_AUTO from SExtractor as the total mag- 
nitude of the detected objects. The SED (colors) was 
measured in the optimal-S/N aperture size for each 
galaxy (see KY05 for details). We excluded regions 
with lower sensitivity, such as the edge region of each 
ISAAC FOV for GOODS-S, from our analysis. The ar- 
eas in the analyses were 94 arcmin 2 for GOODS-S, 3.0 
arcmin 2 for IRAC UDF, 3.4 arcmin 2 for HDF-N, and 
4.1 arcmin 2 for FIRES. We performed simulations us- 
ing the IRAF/ ARTD ATA package to quantify the depth 
of each field. The 90% completeness limit for the point 
source is K ~21.8-22.6 (~21.4 for one FOV) for GOODS- 
S, K ~ 22.5 for IRAC UDF, K ~ 23.0 for HDF-N, and 
K ~ 23.5 for FIRES. 

3. ANALYSIS 

Using these optical and NIR photometric data, we 
measured the redshift, the stellar mass, and the rest- 
frame U — V color of each galaxy in the same manner 
as that in KY05. When possible, we used the spec- 
troscopic redshifts from literature iSzokolv et alJ 12004 
| Le Fevre et"aTI l200l IVanzella et all I 2005L iWirth et all 
12004 ICohen et a.Ul2000l ILa.hhe et a,Ul20ol . For objects 
without spectroscopic identification, we estimated the 
photometric redshifts using the publicly available "hy- 
perz" code (Bolzonella et al. 2000) in the same manner 
as that in KY05; however, for F IRES, we used the pub- 
licly available phot-z catalogue l)Labbe et alJl2003|) . 

After the redshifts were determined, we performed a 
detaile d SED fitting with the GALAXEV synthetic li- 
brary ijBruzual k. Charloti r2003) in order to measure the 
stellar mass and the rest-frame color of the galaxies. 



We assumed IMF of iChabrierl l|2003|L ex ponentially de- 
caying SFR, and the extinction law of iCalzetti et alJ 
(2000). The free parameters were the age, star forma- 
tion timescale, color excess, and metallicity. The fitting 
results were used to estimate the M/L ratio (and subse- 
quently, the stellar mass) and measure the K-correction 
values for the rest-frame colors. The parameter range 
of dust extinction is < E(B - V) < 0.62, and the 
median of the best-fitted values is 0.13 (about 75% of 
the samples had E(_B — V) < 0.3). Additional details 
of the procedures and the discussion on the uncertainty 
of the mass and color evaluation are described in KY05. 
It should be noted that although each parameter such as 
dust extinction, stellar age, or star formation timescale is 
not constrained very tightly, the stellar mass can be esti- 
mated with a relatively low uncertainty. This is because 
the effects of these parameters tend to cancel out in the 
esti mation of the stellar mass, as discussed in KY05 (see 
also lPapovich. Dickinson, fc Fergusonl200il for a detailed 
discussion) . In addition, note that the rest-frame U — V 
color is measured from the observed photometric data 
without extrapolation. Therefore, the stellar mass and 
the U — V color used in the current analysis are robust 
quantities. 

4. RESULTS & DISCUSSION 

In Figure m for each of the four fields, we plotted the 
rest-frame U — V color of the K-selected galaxies in each 
redshift bin as a function of the stellar mass. First, we 
can see from this figure that the color distributions of the 
galaxies in the four fields are very similar. For reference, 
we estimated the weighted mean colors with a bin of 
±0.15 dex of the GOODS-S sample at 0.3 < z < 0.7 (in 
the top- left panel), and plotted them in all the panels by 
dotted-dash lines. The U — V color is always very blue 
for M ste iiar < 10 9 M Q and becomes redder as the stellar 
mass increases toward ~ 10 11 M©. Most of the galaxies 
with M ste iiar > 10 11 M© have the reddest color. 

In fact, it is evident from Figure that over the entire 
redshift range at 0.3 < z < 2.7, the more massive galaxies 
always have the redder U — V color. For example, while 
most of the galaxies with M ste ii a r ^ 3 x 10 9 M© have 
U— V < 0.3 at any redshift, many objects with M ste u ar > 
3 x 10 10 M have U - V ~ I. This trend is consistent 
with the results found in lFeulner et alJ l)2005bl) and other 
previous studies. If we observe the redshift evolution in 
Figure ^ more closely, we notice the gradual bluing of 
the U — V color along the redshift within the mass range 
of 10 9 M < M ste iiar < 10 10 M©. While this systematic 
evolution of the rest-frame color distribution of low-mass 
galaxies has already been reported in KY05 for HDF-N, 
it is now confirmed by the larger sample that includes 
three other different fields. 

In Figure |3 we try to quantify the observed mass- 
dependent color evolution by using an analytical fitting 
function. The data for all the four fields are combined. 
We excluded the lowest mass galaxies in each field from 
the sample so that the color bias near the detection limit 
(solid lines in Figure Q does not affect the fitting sig- 
nificantly; the limiting masses are set so that the K- 
band source detection completeness at U — V ~ 0.8 is 
50%. These limiting stellar masses are tabulated in Table 
2] Considering the shape of the mass-weighted average 
color profile in GOODS-S (the top-left panel in Figure 
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Fig. 2. — Rest-frame U — V color distribution as a function of the stellar mass for the sum of all the four fields in each redshift bin. Data 
from different fields are plotted by different symbols. Blue and red symbols represent the spec-z and phot-z samples, respectively. Solid 
lines show the fitted "tanh" function in each redshift bin. The fitting parameters, Wmin an d Mtr (see text) are also shown. 



TABLE 1 

Lower limit of stellar mass in each field (Mq) 



Redshift 


GOODS-South 


IRAC UDF 


HDF-N 


FIRES 


0.3-0.7 


1. 4-3.3x10* 


1.8x10 s 


1.1x10 s 


6.6xl0 7 


0.7-1.1 


4.2-lOxlO 8 


5.5x10 s 


3.5x10 s 


2.0x10 s 


1.1-1.9 


1.2-2.9X10 9 


1.6xl0 9 


l.OxlO 9 


5.8x10 s 


1.9-2.7 


3.1-7.6X10 9 


4.1xl0 9 


2.6xl0 9 


1.5xl0 9 



1), we chose a hyperbolic tangent as the fittin g function. 
This was in keeping with lBaldrv et al.1 ((2004) who used 
a similar form to quantify the bimodality in the color- 
magnitude distribution of the SDSS sample. The fitting 
function is written as follows: 

U - V(M stella r)=( UVMIN+ 2 UVMAX ) + 
( t/VW-tAW ) x tanh (R(M stellar - Mtr)) 

where UVmin and UVmax are the asymptotic minimum 
and maximum values of the rest-frame U — V color at the 
low- and high-mass ends. M tr and R represent the mass 
and mass range at which the color transition occurs. For 
simplicity, we fixed R and UVmax , which seem to change 
little with the redshift in our sample, to the values fitted 
at 0.3 < z < 0.7 (R = 1.0, UVmax = 1.3). Next, we fit- 
ted the color distribution in each redshift bin by changing 
the values of M tr and UVmin- The result is shown by 
the solid lines in Figure |21 The best-fitted values of the 
parameters are M tr = 1.57^0.23 x 1010 M o and UVmin = 



-0.036 ± 0.055 at 0.3 < z < 0.7, M tr = 1.71±o24 x 10 
M Q and UVmin = -0.277 ± 0.056 at 0.7 <z < 1.1, 
Mtr = 3.56lJ;s4 x 1010 M and UVmin = -0.317±0.128 
at 1.1 < z < 1.9, M fr = 5.77ti ^ x 10 10 M Q and 
UVmin = -0.344 ± 0.084 at 1.9 < z < 2.7. These 
results again show that the U — V color at the low-mass 
end gradually evolves blueward with the redshift. In ad- 
dition, the transition mass also seems to increase along 
the redshift, although t he uncertain t y is re latively large. 
It is noteworthy that iBundv et al.l {2006) also recently 
discovered similar transition mass increments with the 
redshift at 0.4 < z < 1.4 in the study using the DEEP2 
data. 

Figure 01 shows the evolution of the rest- frame U — V 
color at each stellar mass range for the combined data. 
We plotted the weighted average colors binned over the 
redshift range of ±0.15 (dotted-dash lines). In the top 
panel, it is more clearly evident that the color distri- 
bution of the galaxies with M ste ;/ ar < 3 x 10 9 M Q be- 
comes gradually bluer with an increase in the redshift. 
In KY05, we discussed that the star formation activity in 
these low-mass galaxies occurs rather continuously since 
they exhibit a very blue color at any redshift. A bluer 
color even at a higher redshift also suggests younger av- 
erage ages of the galaxies at z =2-3. In the lower panels, 
we see that the colors of galaxies with a larger stellar 
mass also become bluer at a high redshift. There seems 
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Fig. 3. — rest-frame U — V color vs redshift for galaxies within each stellar mass range. The plotted sample and the meaning of the 
symbols arc similar to those in Figure [2] Dotted-dash lines and solid lines show the weighted average and the best-fit GALAXEV models 
with exponentially decaying SFR. Dashed lines show the upper and lower envelope of the color distribution. 



to be a trend where the more massive populations show 
a stronger color evolution at a high redshift, which in 
fact corresponds to the evolution of the transition mass 
found in Figure El 

We try to fit such color evolutions using GALAXEV 
population synthesis models with simple star formation 
histories. We simply consider an exponentially decaying 
SFR and assume solar metallicity and no dust extinction. 
We estimate the U — V color for the model grid with the 
star formation timescale r = 0.01, 0.05, 0.1, 0.5, 1, 1.5, 
2, 3, 5, 15, and oo (constant SFR) Gyr and the forma- 
tion redshift Zf = 3, 4, 5, 6, 7, and 8. Eventually, we 



found that r mainly affects the overall slope of the color 
evolution, and Zf strongly affects the color at a high red- 
shift (e.g., z > 2). The resultant best-fit models and 
their r and Zf values are shown in Figure |3I It is found 
from the results that the more massive populations ap- 
pear to have a shorter star formation timescale. In other 
words, they form the larger fraction of their stellar con- 
tents at a high redshift. Such a trend is very similar to 
the conce pt of mass-dependent star formation histories 
argued bv iHeavens et all ((2004) . Further, the trend is 
also consistent with the studies on mass-dependent evo- 
lutions of specific star formation rates at a high redshift 
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ijFeulner et alJl2005bl Uuneau et all 12005(1 . However, in 
this study, the trend is depicted by more directly observ- 
able quantities. 

In these analyses, it is important to check the effect of 
the color bias near the detection limit. The bias against 
a red high-M/L population is expected to be severe, es- 
pecially for low-mass galaxies at a high redshift since the 
if-band filter samples relatively short rest-frame wave- 
lengths. We finally discovered, however, that our deep 
-fT-selected samples do not suffer considerably from the 
bias. It is true that our 50% completeness limit at 
U — V = 0.8 could neglect some fractions of relatively 
red low-mass galaxies in the highest redshift bin. How- 
ever, the color distribution tends to concentrate at con- 
siderably bluer values (U — V < 0), and there is few 
galaxies that are redder than U — V = 0.3, where the 
completeness is greater than 90%. This can be demon- 
strated by the lowest mass objects near the detection 
limit (K ~ 23.5) in the highest redshift bin of the deep- 
est FIRES data. For example, out of 20 galaxies with 
the stellar mass of 1.5 x 10 9 M Q < M steHar < 3.0 x 10 9 
M Q in the 1.9 < z < 2.7 bin of FIRES, 19 galaxies have 
U — V < 0, and one galaxy has < U — V < 0.1; further, 
no galaxy has 0.1 < U — V < 0.3. 

Indeed, a sample selection with very deep NIR data 
is essential not only for evaluating the stellar mass of 
galaxies but also for generating relatively unbiased sam- 
ples to study the color distribution of these low-mass 
galaxies at a high redshift. For example, /-band selec- 
tions such as those us ed for the FORS Deep Field data in 
iFeulner et alJ l|2005bf) could introduce a se vere color bias. 
The de tection limit of I{AB) = 26.4 in IFeulner et all 
( 2005b) can only sample galaxies with U — V < —0.1 for 
M s t e ;/ ar = 2 x 10 9 Mq at z > 2. Such a selection does 
not reveal whether red galaxies exist. 

In Figure E3 we found that the stronger color evolu- 
tion of massive galaxies at a high redshift was explained 
by the exponentially decaying SFR models with Zf ~4- 
5, for which the star formation timescale decreases with 
mass. On the other hand, the evolution of the average 
color of low-mass (< 3 x 10 9 M ) galaxies can actually 
be explained by the constant SFR model. However, the 
best-fit formation redshift occurs at zt ~ 4, which is sim- 
ilar to those for high-mass galaxies. If this is true, the 
differences in the star formation histories among galax- 
ies with different stellar masses are not mainly in the 
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epochs when the star formation began, but rather, in the 
star formation efficiencies and/or mechanisms by which 
star formation activities are suppressed. The long star- 
formation timescale for low-mass galaxies can be caused 
by the self-regulation o f star formation, for exam ple, by 
a supernova feedback (D ekel fc Birnboiml l2005|) . The 
mass dependence of the star formation efficiencies and 
suppression mechanisms can be related with processes 
such as the effici ent exhaustion of gas reservoirs (e.g., 
iMenci et al.ll2005[) or shock hea ting of cold gas during a 
major merger (|Cox et al.ll2004|) . 

It is noteworthy that some of the assumptions in our 
analysis may be oversimplifications. If we consider dust 
extinction, the estimated formation redshift shifts to 
smaller values. For example, if we assume the extinction 
of Fi(B — V) = 0.13 (median- fitted value of our sample) 
for low-mass galaxies, the best-fit formation redshift of 
these galaxies shifts from z/=4toz/=3. We also 
verified which variations of star formation histories are 
permissible if the scatter of the colors around the aver- 
age values is considered. To illustrate this, we plotted 
the GALAXEV models that trace the upper and lower 
envelopes of the color distribution by dashed lines in each 
panel in Figure |3| 

Finally, we note that several other surveys have investi- 
gated the evolution of the stellar mass or luminosity func- 
tion of field galaxies a nd found that the normalization de- 
creas es at z > 1 fe.g..lFontana et all2004UYamada et alJ 
120051 ICamiti et al .1 1200fl ISa.racco et al J I200rfl . Several 
studies have also su ggested mass-dependent number den- 
sity evolution (e.g.. iGlazebrook et alJl2004L iDrorv et alJ 
12004 iDrorv et alll2005|) . The nextstep is to consistently 
explain both the color (SED) evolution and the number 
density evolution as functions of the stellar mass. 
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